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Based  on  arbitrarily  right  censored  observations,  asymptotically  distri¬ 
bution  free  tests  are  proposed  for  testing  the  equivalence  of  K  treatments 
in  the  analysis  of  a  randomized  complete  block  design  with  random  block 
effect.  Multiple  comparisons  procedures  which  incorporate  a  pairwise  ranking 
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1.  INTRODUCTION 


w"- 


Let  (X^,  .X^)’,  •**  ,  (X^,  ’**  .X^) '  be  independent,  identically 

distributed  random  vectors  with  a  continuous  joint  distribution  function 
and  marginal  distribution  functions  F^  ,  i  *  1,  •••  ,K  .  The  null  hypothesis 
II  ,  which  we  wish  to  test,  is  that  C^  is  symmetric  in  its  K  arguments, 
i.e.  that 

G  (x.  ,  *•*  ,x.  )  =  G  (x  ,  ,x  )  , 

X1  XK  IK 

for  all  permutations  (i  ,  •••  ,i  )  of  (1,  •••  ,K)  .  Note  that  H  implies 

X  K  0 

that 

F°(x)  =  F^(x)  =  ...  =  F^(x)  ,  for  all  x  e  R  . 

Such  a  hypothesis  testing  problem  arises  when  we  test  the  equivalence  of 
K  treatments  in  the  analysis  of  a  randomized  complete  block  design  and 
assume  that  the  blocks  themselves  are  drawn  at  random  from  a  population  of 
blocks.  There  are  several  standard  nonparametric  tests  for  this  case,  for 
example,  Friedman  test  and  the  aligned  rank  order  test  (Mehra  &  Sarangi 
(1967),  and  Sen  (1968)).  For  an  ordered  alternative  : 

F^(x)  <  F^(x)  £  ”•  <  F^(x)  *  For  x  c  R  , 

where  at  least  one  of  the  inequalities  is  strict  for  some  x  ,  we  have  tests 
proposed  by  Doksum  (1967)  and  Hollander  (1967).  If  one  is  not  only  interested 
in  a  single  global  test  of  whether  the  K  treatments  are  equivalent,  then 
multiple  comparisons  procedures  are  more  appropriate  (cf  Miller  (1966,  1977)). 
Now  suppose  that  a  clinical  trial  is  conducted  to  compare  K  treatments. 


the  parameter  of  interest  is  the  length  of  survival.  It  is  common  that  at 


2 


the  end  of  the  trial,  there  may  be  incomplete  survival  information  on  certain 

individuals.  More  specifically,  the  response  X^  may  be  censored  from 

right  by  an  independent  variable  Z  .  so  that  the  n  K-tuples  (X^  ,  •••  ,X?,)' 

ij  ij  Kj 

cannot  always  be  observed.  Rather  our  observations  consist  of  nK  pairs 


^Xij,Aij^'  *  where  Xij 


min(X^,Zij)  and  A^  =  1  ,  if  X^  is  observed 


and  zero  otherwise.  Furthermore,  for  mathematical  convenience  Z, .  is  assumed 

ij 

to  be  governed  by  a  continuous  distribution  function  J^(*)  »  i  ■  1,  •••  K  , 
j  =  1,  •  •  •  ,n  . 

in  the  parametric  statistics,  the  procedure  for  testing  HQ  based  on 
observations  (X^ . ,A  ) ’  is  complicated  and  no  results  have  been  obtained. 
Although  a  generalized  Friedman  test  (Patel  (1975))  is  available  for  testing 
Hq  based  on  censored  observations,  it  is  rather  inefficient  when  there  are  too 
many  censored  K-tuples  in  the  data.  As  an  extreme  example,  for  the  data 
(12+,9+,13) '  ,  (9,8+,5+)'  ,  (4,3+,2+)’  ,  (2,1+,1+)’  ,  where  "+"  denotes  cen¬ 
soring,  the  generalized  Friedman  test  which  uses  information  only  within 


blocks  leads  to  no  conclusion  about  . 

In  this  article,  we  propose  and  analyze  test  procedure!,  which  allow 
interblock  comparisons  for  testing  based  on  observations  (X  ,A^)'  . 

A  multiple  comparisons  procedure  which  selects  treatments,  if  any,  that  differ 
from  one  another  in  terms  of  enhanced  survival  is  also  presented. 

The  corresponding  K  sample  problems  have  been  investigated  by  Breslow 


(1970) ,  Crowley  &  Thomas  (1975),  and  Koziol  &  Reid  (1977). 
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2.  DEFINITIONS  AND  A  BASIC  THEOREM 


—  6  — 

Let  and  F^  be  the  empirical 

tions  of  X. ,  .  Also,  let  F.e  and  F. 

ij  1  i 

subdistribution  functions  of  uncensored 


and  the  theoretical  survival  func- 
be  the  empirical  and  the  theoretical 


F4(b)  =  P(X±j  < s,A  =  1)  and 
~e  1  n 

F.(s)  =  n  X  KX^  -s*An  =  X)  » 
j=l  J  ' 


where  1(A)  is  the  indicator  function  of  the  event  A  ,  i  =  1,  •••  K  . 

Furthermore,  the  joint  subdistribution  functions  of  (X^  ,X^, ^)  '  , 
‘VVj’VVj*'  and  (Xlj ’  **  *  ,XKj)r  are  defined  by 

Dil,(s,t;6,e)  »P(Xlj  <  s,X1,J  <  t *  6,A4, ^  -  e)  , 

4u‘4'4fi'e-4kj'4’>Vj't,)  •  “"d 

®<*1-  •"  ,*Ki{1>  'V  ■  P<K1J  SS1,'"'XKJ  5*K'Alj  "S'  ■"  'Vl’V  • 
respectively, 

where  i  ,  i’  ,  k  ,  and  k’  denote  distinct  indices.  The  corresponding 
empirical  subdistribution  functions  are  defined  by 

(s,t;6,e)  =  n"  KX^  5  s.X^  5  t.A^  *  6,6^  -  g)  , 


Hii,kk'(s,t’S,’t,;6’e,,S'’e’) 


-1 


£ 1<xu s  s-xi'j  j  •  ' 4 ’-V3  ■  E'> 


ir  Vi  I 
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and 

~e 


(Sl*  **’  ,SK;61*  **  ’  ,5K>  “  n  ^  I(Xlj  ~  ®1»  "*  ,XKj  "  SK’ 


At4*6i»  *•*  -6_> 


1J  1 


Kj  K 


All  the  procedures  proposed  in  this  article  are  based  on  the  Cehan  scoring 
function  ip  ,  where,  for  any  two  pairs  of  observations  (X  ,A^  ) '  and 

(X^ i j i , A^ i j | )  , 


^<Xij,AlJ*Xi,j,»Ai,j,) 


/  +  1  ,  Xy  >  XlfJ,  ,  .  ~  1 

-1  •  Xij  *  Aij  =1 

0  ,  otherwise. 


Now,  let  pu,  *  E^<Xik’Aik’Xi'  k'  ’^i'k*  >  ’  where  1  *  and  k  t  k'  • 

Under  HQ  ,  p^,  *  0  which  is  independent  of  censoring  distributions  and 

J^,  (cf  Efron  (1967)  and  Mantel  (1967)).  Consider  Gehan's  statistic  (1965) 

V,,.  for  populations  i  and  i’  : 
ii,n 


_3  n  n 


vii;n  •  "  *  ^  'HWVr'\T) 


The  following  Theorem  whose  proof  is  in  the  Appendix  provides  the  asymptotic 
distribution  of  random  vectors  <V*  >  ,  where  V* 


-4  n  n 


ii,n 


ii,n 


"  j,E=1('('(Xij,Aij’Xi,j,,Ai'j,)"Pii’)  ’  i’1'  "  lf  *“  ,K  80d  1  *  *’ 


Theorem  1.  As  n  -*•  ®  ,  the  random  vector  <vli»n>  converges  in  distribution 
to  a  multivariate  normal  random  vector  with  mean  0  and  covariance  matrix 

rs* 

t  *  (ollf  kk,)  ,  i  <  i'  ,  k  <  k'  .  Under  Hq  , 


°ii’  ,kk*  "  aii\kk' 
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E<V(  V  -  Vl'(  V  +  V3Fl(Xl\|>  -  VVi>) 

(?k'‘*kj>-4kj?ki(xkj)+VA(xk-j,-VVj))  • 


If  the  censoring  variables  ,  ***  ,Z^  are  interchangeable  (hereafter 

referred  to  as  Condition  A)  ,  then  J,  =  J.  =  •••  =  J,,  and,  under  H  ,  the 

1  L  K-  U 

random  pairs  (X  ,A  ) *  •••  ,  (X..,,A  )'  are  also  interchangeable,  j  =  1,  •••  ,n 

lj  iJ  .  KJ 

For  future  references,  let  F  =  F,  =  *•*  =  and  F  =  F  =  •••  =  F„  and  also 


let 


o2  =  e(f(X  .)  -  A. ,F(X  .))2  and 

\  ij  lj  ij  / 


t  -  E(?(*1J)-4yF(Xy>)(F(X1.j)-41,jr<X1,j)) 


3.  MULTIPLE  COMPARISONS  UNDER  HQ 

When  it  is  less  important  to  test  the  overall  hypothesis  than  to 

compare  certain  populations  on  an  individual  basis,  a  multiple  comparisons 

procedure  is  more  desirable  than  a  single  global  test.  The  decision  of 

which,  if  any,  populations  differ  from  one  another  can  be  made  by  examining 

the  individual  comparison  V..,  .  We  would  conclude  that  populations  i 

ii,n 

and  i'  differ  if  lv^j»n!  *s  sufficiently  large.  A  conservative  multiple 
comparisons  procedure  at  level  a  can  be  obtained  through  an  inequality  due 
to  Si dak  (1967),  i.e.  if  U  =  (U^,  ,Ur) '  is  normally  distributed  with 
mean  J)  and  correlation  matrix  R  ,  then 
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P  (  max  ju J  <  u  |  U~N  (0,R))  * 

l<i<r 

P  (  max  | U  |  £  u  |  U  ~ N  (0,1  ))  . 
r  l<i<r  1  r  ~  r 


It  follows  from  Theorem  1  that,  for  large  n  ,  one  would  judge  that  populations 
i  and  i'  as  different  at  level  a  if 

'vii,,J/(aii',ii’  )2  >  ?B/2  ’ 


where  is  t*le  100(1 -P/2)  percentage  point  of  the  N(0,1)  distribu- 

K (K-l) /2  *0 

tion,  (1-6)  =  1-a  ,  and  a..,  is  a  (strongly)  consistent 

ii  ) xi 

estimator  of  : 


1  1 


ii’ , ii' 


=  22  f  [F*,(s) -6F*,(s) +EF?(t) -F*(t)]2dD?  ,(s,t;6,e) 

6=0  e- 0  J  1 


Another  multiple  comparisons  procedure  can  be  obtained  through  the 

Bonferroni  inequality.  But  as  indicated  by  Koziol  and  Reid  (1977,  p.  1154), 

this  procedure  is  less  powerful  than  the  one  based  on  &idak  inequality. 

Some  improvement  of  the  above  results  can  be  made  under  Condition  A  (the 

interchangeability  of  Z  ,  •  *  *  ,ZV.)  .  It  can  be  shown  that,  under  H  and 

lj  y  ,  u 

Condition  A,  the  asymptotic  covariance  structure  of  < — 4^  *!;-  >  is  identical 

(az-r)*S 

to  that  of  the  vector  having  K(K-l)/2  components  Y^-Y^,  ,  1  s  i  <  i’  ^  K  , 

where  Y  *  (Y, ,  •••  ,Y„)  ’  ~  N,.(0,Tv)  ,  and  I  is  the  identity  matrix.  There- 
~  1  l\  K  ~  K  K 

fore,  one  would  conclude  that  populations  i  and  i’  differ  at  level  a  if 


I’iC  J  ,  a 

(a2-f)f  K 
n  n 
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where  q  is  the  lOO(l-ct)  percentage  point  of  the  distribution  of  the  range 
K 

of  K  independent  unit  normal  random  variables  (cf.  Pearson  and  Hartley  (1966), 

~  2 
Tables  22  and  23)  and  o  and  T  are  (strongly)  consistent  estimators  of  o 

n  n 

and  t  : 

o2  =  /  "  Fe(s))2dFC(s)  -  J( FC(s))2d(F(s) +  FC(s))  , 

/v\  1  K-l  K 

(,)t  =  S  S  S  /(Fe(s)  -6Fc(s))(Fe(t) -eFe(t))dD®.,(s,t;6,e)  , 

V2/  n  6,e=0  i-1  i’=i+l  J  11 

Fe(s)  *  K"1  y  Fe(s)  and  F*(s)  =  K_1  £  F^s)  . 


4.  TESTING  HQ  AGAINST  ORDERED  ALTERNATIVE 

In  this  section,  we  present  a  test  statistic  which  is  a  generaliza¬ 
tion  of  the  Jonckheere  K-sample  test  (Jonckheere  (1954),  Patel  &  Hoel  (1973)) 
for  testing  against  based  on  observations  (X^,A^)'  ,  where 


S 

n 


K-l  K 

l  ,1 

i*l  i  =i+l 


iir,n 


Large  values  of  suggest  rejection  of  the  null  hypothesis  .  It 

follows  from  Theorem  1  that,  under  H„  ,  S  converges  in  distribution  to  a 

On 

normal  random  variable  with  mean  0  and  variance 

K-l  K  K-l  K 


■'  A  |\  « 

2  2  2  2  kk. 
i«l  i'-i+l  k-l  k'-k+l  1  ,K 


Note  that  o^,  *  where  i  ,  i’  ,  k  and  k*  are  distinct  integers,  can  be 

estimated  consistently  by 


8 


-0 

ii’.kk' 


I  /"(f®,  (s)  -<SF*f(s)  +  e??(t)-F?(t)) 

l,e,6’ ,e'=0  J  V  1  1  1  17 


(f®,(s')  -6,F^,(s,)  +  e,F^(t') -F^(t'))d  Mj1*kk*(s.t.s,,t'  ;6,e,6\e’) 


5.  A  GLOBAL  TEST  FOR  11Q 

Tf  we  are  interested  in  testing  HQ  against  a  general  alternative 
hypothesis  that  G°  is  not  symmetric  in  its  K  arguments,  then  a  global  test 
is  more  appropriate.  Unfortunately,  it  is  not  feasible  to  get  such  a  test  under 
unequal  patterns  of  censorship  operative  on  K  populations.  In  this  Section, 
we  assume  that  the  censoring  variables  ,  •  •  •  .Z^  are  interchangeable 
(Condition  A)  and  present  a  generalized  Breslow  test  (1970)  for  testing  HQ 
based  on  observations  (X^.A^)’  . 

Define  a  vector  score  statistic  W  =  (W  ,  •  •  •  W  ) '  ,  where  W  is  the 

~n  In  Kn  in 

total  score  comparing  the  ith  sample  with  the  remaining  (K-l)  samples,  i .e. 

K 

w,  -  y  v.., 

in  ii'.n 

i’*i 


It  follows  from  Theorem  1  that,  under  H_.  ,  W  converges  in  distribution  to 

u  ~n 

a  normal  random  vector  with  mean  0  and  covariance  matrix 


where  y 


n  -  , 


r  - 


A 
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4> 


_  K  _ 

(K-1)(F(X1.)  -  &1.F(X1.))  +  V  A  .F(X  ,)  and 
lj  lj 


♦'  •  0=-l)(«X1.j)-4KjF<XKj))  +  |y|ljF<X|lj)  . 


2  ^2 
Again,  y  and  n  can  be  estimated  consistently,  for  example,  by  y^  and 

A 

n  ,  where 
n 


and 


n 

n 


s  /r<K-ix»*(.1)-«1?(.1))+ 2  «.?(..)] 


m=2 


K-l 


[(K-l)(Fe(sK)-5KFe(sK)  +  y  6^(8^)]  dGe(s1,---,sK;61,---,'SK) 

£=1 


2 

The  latent  roots  of  T  are  (y  +  (K-l)h)  ,  which  is  zero,  with  the  latent 

vector  q  =  (1, • * • ,1) ' /K^  and  (y^-n)  of  multiplicity  (K-l)  with  (K-l) 

latent  (any  set  of)  vectors  orthogonal  to  (1, •••,!)'  »q^*’'",q^  ^  ,  say 

(cf  Rao  (1973),  p.  67).  For  example,  when  K  =  3  ,  =  (y)^  (-^,1 ,-%) '  , 

q2  *  (j) r(%,-%»0) '  ,  q3  =  (— )*(1,1,1)'  .  In  general,  the  q±  are  perhaps  most 

easily  found  by  means  of  the  Gram-Schmidt  orthogonalization  process.  Since  T 

is  symmetric,  therefore,  Q'TQ  **  ,  where  Q  =  (q  ,  •••  ,q  )  which  is  an 

2 

orthogonal  matrix,  A  -  diag  (A,  •••  ,A,0)  and  A  ■  y  -  f|  .  It  follows 

K  x  K 

that,  under  , 


T 

n 


.  K-i 

A  V  (qjw  r 
n  ^  v*i~n 


is  asymptotically  distributed  in  a  chi-squared  distribution  with  K-l  degrees 
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A  S\  2  A 

of  freedom,  where  A  =  y  -  n  •  Large  values  of  T  suggest  a  rejection 

n  n  n  n 

decision  of  Hq  . 

For  K  =  2  and  unequal  patterns  of  censorship,  an  asymptotically  distri¬ 
bution  free  test  was  studied  by  Wei  (1979)  for  testing  based  on  observa¬ 

tions  (X  ,A  )'  ,  i  =  1,2  ;  j  =  1,  ,n  .  This  test  is  reduced  to  the 
above  generalized  Breslow  test  when  the  two  censoring  distributions  and 

are  equal. 


6 .  REMARKS 

In  practice,  all  the  consistent  estimators  in  this  article  can  be 
obtained  through  simple  summation  processes.  Therefore,  it  is  feasible  to 
evaluate  the  test  procedures  by  hand  calculations  and  extremely  easy  if  a 
small  programmable  calculator  is  available. 

Although  there  are  several  multivariate  parametric  models  available  for 
life  testing  in  reliability  (cf  Freund  (1961),  Marshall  &  Olkin  (1967),  and 
Block  &  Basu  (1974)),  none  is  satisfactory  for  survivability  theory,  espe¬ 
cially  for  our  current  settings.  For  example,  there  is  no  continuous  multi¬ 
variate  distribution  function  which  has  exponential  marginals  and  possesses  a 
physical  interpretation  in  survival  analysis.  This  probably  makes  nonpara- 
metric  procedures  more  appealing  in  multivariate  survival  studies. 

Also,  for  lack  of  appropriate  parametric  models,  the  efficiency  study  of 
our  new  test  procedures  seems  unfeasible.  However,  since  the  new  procedures 
allow  interblock  comparisons,  we  expect  that  they  are  more  powerful  than  the 
procedures  based  on  Friedman  statistic. 
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7.  APPENDIX 


Proof  of  Theorem  1. 


=  E^(Xij,Aij’Xi'j,,Ai,j,)lXi,j',Ai'j']  "  Pii'  3nd 
*ii,(Xl,j,»Al,j,)  =  Ef*<Xij  *Aij  *X1*  j  *  *A1' j  *  lX£J  *Aij  1  "  Pil'  '  Where 


Also,  let 


V**, 
ll  ,n 


n-%(j!(<1,<Xlj,41.)+<.;i,(Xi,J,A1,J))) 


8(Xij  ,Alj*Xl' j'  »Airj»)  '^Xij  ,Alj  *Xi*  j  *  ,Ai’  j  *  ^  ~  Pii'  ~ 
^ii’ (Xij  *Aij)  ~  ^il* (X1*  j ' ,A1' j  » }  * 


Then, 


E(V*  ,  -V**  )' 

xi  ,n  ii  ,n 


-'Ml,  X  g<vvxi,j,’Vj,)] 


n  n  n  n 


=  n"  I  I  1  1  h(j  ,j  ’  , 

j=l  j’-l  £.=1  £'=1 


(A.  1) 


where  h(j,j  ,£,£')  -  E  [K  ( j  j  i  j  '^8^Xi£,’^i£’Xi»£i  *Aj_  i  ^>)  1 

Because  g  is  bounded,  we  can  ignore  any  u  terms  of  the  sum  in  (A.l)  if  u 

3 

is  of  order  0(n  )  .  Consider  the  following  cases  for  which  the  number  of 

3 

terms  is  with  order  larger  than  or  equal  to  0(n  )  (where  j , j  ,£,  and  £,' 
represent  four  distinct  indices) : 


12 


(1)  h(j»j  )  ~  E  g(X^  ’X^f  j  » »A^t  j  '  )E  8(X^  ,AiJ!,,Xi 1  £ » »A^  *  £*  ^  a  ®  * 

(2)  h(j ,  j  ,£,£')  -  h(j,j',*,JD  =  0  ; 

(3)  h(j,j',j,r>  =  h(j,j',£,j») 


=  E{E[g(X±J  »Aij.X.,  ^  ,  ,A1,  j  ,)g(Xlj  »Ai:j  .xil£.  lXij  «Aij  ^ 

=  E  J  E  [g(  X,  j  »A±j  *  X± » j  »  j  ’>  I  X-jj  »A1j]^rfi(Xij  *A^j  *Xl*£>  *^1'  S,f  I  X±j  »A^j]  I 


(A)  h(J,J'  ,*,j)  =  Mj.j’.j’,*')  =  EE[g(Xlj,Aij,Xi,jl,Ai,j()g(X1£,Au,Xi,  ,Ai, 

°  ElE[g(X.j,Aij,Xi,j,,Ai,.,)|xij,Aij,Xi,j,Ailj] 

E[8(XU’Au»Xi»j»Ai*j)lXij»Alj.Xi.J.Ai,j]}  =0  . 

rx 

It  follows  that  E(V*  ,  -V**,  )  -+  0  ,  as  n  00  .  By  Corollary  6  of 

xi  i n  n  jit 

Lehmann  (1975),  p.  289)  and  the  fact  that  the  vector  <V**,  >  has  a  multi- 

ii  ,n 

variate  normal  limiting  distribution  N(0,t)  ,  the  vector  < V*  .  >  has  the 

~  ii  ,n 

same  multivariate  normal  limiting  distribution  N(0,$)  .  Under  ,  p^,  = 

it  follows  that  V .  . ,  =  V*  ,  , 

ii  ,n  ii  ,n 


V°W  •  W  -Vi'(V  and 


^iir  j 1  j  f )  *  j  '  j 9 )  ~  * 


Q.E.D. 
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